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Jet charge has played an important role in experimentally testing the Parton Model and the 
Standard Model, and has many potential LHC applications. The energy-weighted charge of a jet 
is not an infrared-safe quantity, so hadronization must be taken into account. Here we develop 
the formalism to calculate it, cleanly separating the nonperturbative from the perturbative physics, 
which we compute at one-loop order. We first study the average and width of the jet charge 
distribution, for which the nonperturbative input is related to (dihadron) fragmentation functions. 
In an alternative and novel approach, we consider the full nonperturbative jet charge distribution 
and calculate its evolution and jet algorithm corrections, which has a natural Monte Carlo-style 
implementation. Our numerical results are compared to Pythia and show reasonable agreement. 
This calculation can directly be extended to similar track-based observables, such as the total track 
momentum generated by an energetic parton. 



I. INTRODUCTION 

Jet charge has played an important role in studying 
and testing aspects of the Parton Model and of the Stan- 
dard Model. It was first proposed as a way of measuring 
the charge of a quark Q] . To reduce the sensitivity to ex- 
perimental noise, a weighted-definition of jet charge was 
proposed in Ref. [2], 

Q« = E z ZQ h > (i) 

where the parton i initiates the jet, and Zh = E^/E is the 
fraction of the jet energy carried by the hadron h with 
charge Qh- Various choices for k between 0.2 and 1 have 
been considered at experiments, and we will comment on 
the optimal choice in this paper. Ref. [5] also presented a 
calculation of the jet charge in a recursive, probabilistic 
model. 

A first experimental application of jet charge was in 
deeply-inelastic scattering (DIS) experiments [3H3, find- 
ing evidence for quarks in the nucleon. Jet charge 
also played a role in the measurement of the forward- 
backward asymmetry in e + e~ collisions [H[7], which tests 
the electroweak sector of the Standard Model. Here jet 
charge was used (roughly speaking) to distinguish the 
quark and anti-quark jet. Other applications of jet charge 
include: distinguishing b from b in neutral B meson os- 
cillation measurements [8], assigning jets to the hadron- 
ically decaying W + W~ in a measurement of the triple- 
gauge-boson couplings [9], and excluding an exotic top 
quark model at the Tevatron [TU] . 

We anticipate that jet charge will play a role at the 
LHC, since it is almost the only handle we have to sep- 
arate quark from anti-quark jets or (more ambitiously) 
distinguish quark flavors [11J . Since jet charge is a track- 
based observable it is relatively clean from e.g. pile-up 
contamination. Calculations of jet charge have so far re- 
lied on Monte Carlo programs such as Pythia (T2l [13] . 
Here we calculate the jet charge using and extending 
the framework for fragmentation in jets [TIhTU] in Soft- 



Collinear Effective Theory (SCET) [20H23]. A brief 
overview of some of our main results was presented in 
Ref. [TT]. An advantage of our approach is a clean sep- 
aration of the perturbative physics, described by sim- 
ple analytic formulae, and the nonperturbative physics, 
which needs to be modeled or extracted from data. In 
addition, it provides an estimate of the uncertainty and 
is systematically improvable by calculating higher orders 
in perturbation theory or including power corrections. 

We want to point out that jet charge is not infrared 
(IR) safe, as the following example illustrates. A jet 
consisting of a single quark has jet charge Q q K = Q q . 
If it radiates a gluon with energy fraction z, it has jet 
charge Q q K = (1 — z) K Q q . However, in the collinear limit 
these two configurations arc indistinguishable. Since 
Q q 7^ (1 — z) K Q q , this is not IR safe (unless n = 0). Jet 
charge must therefore be defined at the level of hadrons, 
and hadronization effects are crucial. This is also ob- 
served in Pythia, where the hadronization corrections 
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FIG. 1. The jet charge distribution in Pythia at the hadronic 
and partonic level, for a it-quark jet with E = 100 GeV, using 
the e + e~ anti-fcr algorithm with R=0.5. 
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to jet charge are large, as shown in Fig. [T] 

In this paper we follow two approaches. First we cal- 
culate the average and width of the jet charge distri- 
bution. These arc the two most important quantities 
that determine the power with which jet charge can be 
used experimentally to separate quarks from anti-quarks 
and distinguish quark flavors. Fig. [T] illustrates the chal- 
lenge: the average jet charge is small compared to the 
width. We find that the average charge of a quark jet 
is given by a perturbative coefficient, that contains the 
dependence on the jet algorithm, multiplied by a nonper- 
turbativc number, that can be related to fragmentation 
functions (FFs) [2H US]- The situation is similar, but 
more complicated, for the width of the jet charge dis- 
tribution. The width of quark and gluon jets mix and 
dihadron fragmentation functions now also contribute to 
the nonperturbative input. 

To go beyond the average and width of the jet charge, 
we pursue a second approach where we start with a non- 
perturbative charge distribution. We determine the evo- 
lution and the jet algorithm corrections to this distribu- 
tion at one-loop order. On taking the appropriate mo- 
ments, this reduces to our first approach. A novel feature 
of this second approach is that the equations can be nat- 
urally solved through a parton shower plus hadronization 
model. 

We will show numerical results for the average and 
width of the jet charge, and compare with Pythia. This 
comparison will be both at the nonperturbative level, 
where we take FFs as input, and the perturbative level, 
where we take Pythia as input and compare the calcu- 
lable dependence on the jet energy E and size R. The 
uncertainties on charge-separated FFs are fairly large, so 
the former comparison is not particularly constraining. 
We find reasonable agreement, suggesting that Pythia 
suffices for initial studies of jet charge. Of course this 
may change as the precision of FFs or the knowledge of 
jet charge progresses. 

There are other track-based observables to which our 
framework can be applied, such as the number of charged 
hadrons (tracks) in a jet. Here one could again impose 
a z K -weighting as in Eq. ([!]), or alternatively a cut on 
z, to remove soft radiation and obtain an experimentally 
viable quantity. For k — 1 this corresponds to the total 
track momentum generated by an energetic parton. 

In Sec. |TIJ we calculate the average and width of the 
jet charge distribution, and discuss the relationship with 
(dihadron) FFs. We also introduce dihadron fragmenting 
jet functions here, and discuss some of their properties. 
The approach involving a nonperturbative jet charge dis- 



tribution is presented in Sec. Ill and its Monte Carlo im- 
plementation is described in Sec. |IV| Numerical results 
and a comparison with Pythia are contained in Sec. |V] 
Here we also discuss the optimal choice for k. In App. |A| 
and App. [EJ we give the perturbative and nonperturba- 
tive coefficients that are relevant for our calculation. 



II. AVERAGE AND WIDTH OF THE 
JET CHARGE DISTRIBUTION 

A. Average Jet Charge 

We start out by calculating the average charge for a 
quark jet. (A gluon jet always has average charge zero.) 
The average of the jet charge in Eq. ([!]) is 



Cg-jet J 

!\ zz «y,qh 1 d(Tfe r jet , (2) 

JQ h °q- 



1 -jet 



where cr g _j 0t is the cross section for producing a quark 
jet, and aheq-jct is the cross section for producing the jet 
in which a hadron h is observed. 

We calculate these cross sections in the framework for 
jet production of Refs. . For simplicity we only 

consider e + e~ collisions and use the e + e~ version of fcy- 
type algorithms, with jet size R. (For narrow jets, R <C 1, 
the extension to pp collisions simply amounts to replacing 
the jet energy E by the jet transverse momentum pt, as 
discussed in App. A 2 ) A jet energy veto A is imposed, 
to restrict the radiation in the region between jets. The 
framework is valid for sufficiently narrow, well-separated, 
energetic jets 



tan 2 (i?/2), 



A 



tan 2 (i?/2) 
tan 2 (V>/2)' En 



< 1, 



(3) 



where i/j is the minimum angular separation between jets 
and E m - m the minimum jet energy. It is worth empha- 
sizing that the assumptions in Eq. ([3]), which are used 
to derive Eq. @, are more restrictive than is really nec- 
essary. For example, our results will still hold when two 
jets become closq^J as long as these are not the jets whose 
charge we want to determine. 

Using factorization in SCET, the cross sections at lead- 
ing power are (schematically) given bjj^] 



N-l 



d$ N tr[H N S N ]( J] J^Ji{E,R,ij), (4) 
<l*he fl -iet_ l d ^ Ntr[HNSN] (^ Y[ J t )g^E,R lZ ,n). 



The massless TV-body phase-space for the jets is denoted 
by d<I>Ar. The hard function Hn describes the hard pro- 
cess, and the soft function Sn the soft radiation. Both 



1 The modifications of Eq. |4jl for this case can be obtained from 
Ref. |29j . and do not affect our results. 

2 This involves the factorization of the phase-space constraints 
from the jet algorithm into separate restrictions on the soft and 
collinear radiation, as discussed in e.g. Refs. [27 30 . 
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Hn and Sn are matrices in color space, and the trace is 
over color. For each of the jets there is a jet function Ji 
describing the collinear radiation in the jet. To simplify 
the discussion, we have singled out a jet of flavor i with 
energy E. 

When a hadron is observed in this jet, Ji is replaced 
by a fragmenting jet function (FJF) Q% [3 [T51 [TS] . We 
will neglect the contribution from soft radiation to the jet 
charge, which is suppressed by 0(X 2k ), where the size of 
the SCET power counting parameter A is set by Eq. ^ . 
Since soft gluons do not produce an average charge, this 
additional approximation is absent for the average jet 
charge. We briefly remind the reader of the most impor- 
tant properties of the FJF: has the same renormal- 
ization as the jet function Ji, and can be matched onto 
fragmentation functions D^, 

G?(E, R, z,v) = J2 [ ~ME, R, z', n)D*(^,v) 

J z z \ z / 



X 1 + 



A 2 



AE 2 tan 2 (i?/2) 



(5) 



The perturbative Jij contain the jet algorithm depen- 
dence. The nonperturbative FF describes the frag- 
mentation of an energetic hadron h from a parton j in 
inclusive processes (i.e. without a jet restriction) [24]I25|. 
The F JFs also satisfy certain sum-rules [131 HE] • 

Inserting Eq. Q in Eq. we find that most pieces 
cancel in the ratio of cross sections, yielding [TT] 



(OS) 



g*(E,R,z,n) 
2(2n) 3 J q (E, R, 



JggjE, R, K, fl) 
2(2n)3J q (E,R,ri 



Y^QhD^k^). (6) 



Thus jet charge is independent of the process, up to the 
power corrections in Eq. ([3]). This also implies that jet 
charge is not sensitive to non-global logarithms [3TJ [33] 



in the soft function [33 35 . Because the FJF and jet 
function have the same anomalous dimension, we see that 
the /i-dependence cancels in Eq. Q, as should be the 
case. The last line was obtained using the matching in 
Eq. ([5| for the /t-th moment, with 

Jij{E,R,K,n) = f dzz K Jij(E,R,z,(J,), 
Jo 

ty{K, t i)= [ Azz K D h q {z^). (7) 

As Eq. ^ shows, the average jet charge depends on one 
nonperturbative number (for each quark flavor and k) 



(8) 



related to FFs. Values for at /i = 1 GeV are given in 
App.[B]for several FF sets. 

The matching coefficients J7y and jet functions Ji that 
enter in Eq. ^ are given at next-to-leading order (NLO) 
in App. [A] for /c^-type algorithms. (Results for cone al- 
gorithms can be found in Ref. [H].) Since there is no 
distinction between the various fcT-like jet algorithms at 
this order, we study the difference between the average 
and width of the jet charge distribution in the e + e~ ver- 
sion of hr [55] . Cambridge- Aachen (C-A) [371 I2H] and 
anti-fcr [33] in Pythia. The results are shown in Fig. [2] 
As these jet algorithms only differ in how they cluster 
soft radiation, it is not surprising that the difference is 
only a few % and grows for small R. The dependence 
on k is counterintuitive: one would expect better agree- 
ment for larger values of k, since that suppresses the soft 
radiation. Of course this effect is rather small. 

The matching coefficients J7ij and the jet functions Ji 
contain logarithms of 2Etaxi(R/2)/[i ) so one should take 
fi ~ 2Etaxi(R/2) ~ ER to avoid large logarithms. That 
this combination of E and R appears can be seen by 
boosting the jet along the jet axis. This boost invari- 
ance is of course spoiled by soft radiation from other 
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u-quark, anti-Jfcr, k=1, 2E tan(R/2)=100 GeV 
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FIG. 3. The variation of the average and width of the jet 
charge distribution, keeping 2_Etan(i?/2) = 100 GeV fixed. 



jets, which is accounted for by the power corrections in 
Eq. |[3]). We have investigated to what extent Pythia 
shows the same invariance in Fig. [3j which provides an 
estimate of these power corrections. As you can see, the 
invariance holds at the percent levej^j 

To evolve D®{n,n) to n ~ 2£tan(i?/2), we need its 
RGE 

M) = ^P w («)5«(«, /i) , (9) 

which follows directly from that of the FF. Explicit ex- 
pressions for the one-loop splitting functions are given in 
App. [XJ Note that the mixing with gluon FFs vanishes 
to all orders, since a fragmenting gluon produces no net 
charge. At leading order, the perturbative coefficient in 
Eq. (6| is 1 and the E and R dependence of the jet is 
simply governed by Eq. (p)]) . This shows that the average 
jet charge reduces (dilutes) for larger values of E and R. 
Of course the relevant question is how this changes rela- 
tive to the width of the jet charge distribution, which is 
what we turn to next. 



B. Width of the Jet Charge Distribution 

The width of the jet charge distribution is given by 

(K) 2 = {(Qi) 2 )-{Qi) 2 . (io) 

It is the same for quarks and anti-quarks and does not 
vanish for gluon jets. In the previous section we calcu- 



3 This is process dependent and is expected to be larger for 
hadronic collisions or when more jets are present. We have also 
taken a quite strong cut A on radiation outside the jets here. 



lated (Q l K ), so now we will determine 

(Wk) 2 ) = E E s hu ... >hn fd Zl ---dz n (ii) 

ml*. f, J 



n h\ : ...,h n 



2 1 d"o7 tl .../ tre g;_j et 
(Ji- jet dz\---dz n 



where Sh x ,...,h n is the appropriate symmetry factor. After 
integrating over most of the Zj this simplifies to 



1 do-/, 



'i-jet 



i£z— jet 

dz~ 



Jo h a 

\ dzi dz 2 z^z^y2Q hl Q h 

JO Jo , u . 



+ 



hi,h 2 
1 dcr/ lli / l2 gj_j e t 

o-j-jot dzi dz 2 



(12) 



Here we used Shh = 1/2 and accounted for the factor of 
2 for h\ h 2 by including both orderings. The calcula- 
tion of the first t erm is completely analogous to that of 
(Q l K ) in Sec. II A For the second term we get a dihadron 
FJF g£ lh2 instead of a single-hadron FJF. The relevant 
properties of the dihadron FJF (discussed in the next 
section) are that it also has the same renormalization as 
the jet function Ji , and that it can be matched onto (di- 
hadron) fragmentation functions according to Eq. (jTsj) . 
Using this, we find 

<(Q«) 2 ) 



2(2 7 r) 3 Ji(£,U, / u) 

f dzi f dz 2 z^z% V Q hl Q 
Jo Jo ,,. , 



G^(E,R,z 1 ,z 2 ,») 
2(2tt) 3 Ji(E,R,fi) 



E 



J l3 {E,R, 2k, n) D hlh2 (n k u) + 

2(2%)* J^E, R,fi) L^ Qh ^ h * D i ^ K ^) + 



hi ,h-2 



Jij(E,R,K,fi) 
2{2ir) 3 Ji(E,R,n) 

= E 



E QhiDj 1 (k, v)Qh 2 D'* 2 m (n, n) 



J l3 {E, Rl 2n^) r 5 Q 2 ( 2K , M )+5^( K , M )] 



2{2vf J % {E,R^) 



J gq (E,R,K,iJ,) V^rfiQ, ,,)]'' 



H ' B {2-KfJ g {E,R,^) 



(13) 



as reported in Ref. |TT] (except for the last term that only 
contributes for gluon jets). Here 

J ij (E,R,K,fj l )= [ dzz K (l- zYJij{E,R,z, /x), (14) 
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and the flavor index a(ij), used in an intermediate step, 
is given by 



a{ij) = { q if ij 
9 



qg,gq, 
qg,gq, 
99, qq, qq, 



(15) 



where the q's and g's are of the same flavoiQ The last 
line in Eq. ( 13 ) only contributes for i = g, and leads to 



nonlincarities. These arc not important due to the small- 
ncss of J gq (see Table [v}, except at very small values of 
\x or k. 



Eq. (13) involves two new nonpcrturbativc numbers 



Df (2k, M ) = ^Q^(2k, M ), 

h 

Df Q (K,fi) = Qh l Qh,Df h2 {^^y) 1 (16) 

hi ,h,2 

where the latter depends on dihadron FFs. To evolve 
these to the appropriate scale ~ 2i?tan(i?/2), we use 



fi—D 



QQ 



(Xm) = E 



(17) 



E 



Mm) 

7T 

a s (n) 



which can be obtained by taking moments of the di- 
hadron FF evolution in Eq. (23). The combination 

Df (2k, jit) + Df Q (n, /i), which appears in Eq. ( 13 ), sat- 



isfies this same RG equation. The last line of Eq. (17) 



only contributes for i = g, as in Eq. (13). This nonlinear 



contribution to the evolution is again negligible, unless /i 
or k is very small. However, the mixing between quarks 
and gluons is quite sizable, as we will see in e.g. Fig. [Sj 



The dihadron FJFs can be matched onto (dihadron) 
FFs pi] 



2 ^ R Z\,Z2,H) 



(18) 



J u z ■> \ u u 



E 

j-k 



du dv 



u v 



J ljk (E, R, u, v, n)D^ m)^ 2 , ft) 



In the first term the hadrons fragment from the same 
parton j, and the coefficients Ji 3 are the same as in 
Eq. (|5j. The nonperturbative dihadron FF D^ lh ' 2 de- 
scribes the fragmentation of energetic hadrons hi and h% 
from a parton j in inclusive processes (without a jet re- 
striction) [30HI2] . In the second term of Eq. (fl~8|) the 



hadrons fragment from different partons j and k, which 
is described by matching onto two FFs, and does not con- 
tribute at tree- level, J^l — 0. Performing the matching 
by replacing the hadrons by partons, we find at one-loop 
order that 

J$(E, R, u, v, fj,) = Ji ( p{E, R, u, n)S(l - « - v)6 kMij) . 

(19) 



The 5k, a (ij) indicates that at one-loop order the flavor k is 
completely fixed by ij, as described by a(ij) in Eq. (15). 
The dihadron FJFs satisfy the sum rule 



E [dz 2 z 2 g^ h2 (E,R,z 1 ,z 2 ,v) 



= (1 - zi)g i 1 (E,R,zi,z 2 ,(J,) 



(20) 



which follows from momentum conservation. 

We now perform some basic checks. First we note that 
by using Eq. ( |20[ ) and the corresponding sum rule for 
dihadron FFs |40ll42] 

£ Jdz 2 z 2 D> 3 l ^(z 1 ,z 2 , f i) = (l-z 1 )D^(z 1 ,n), (21) 



C. Dihadron Fragmenting Jet Functions 

The dihadron FJFs g^ lll2 (E,R,z 1 ,z 2 ,fi) are a direct 
extension of the single- hadron FJFs of Refs. [TH [T51 UH] , 
where an additional hadron is observed in the jet. This 
modification of the (IR) state, does not affect the (UV) 
renormalization, which is therefore the same as the jet 
function Ji. 



4 This is true at one loop and requires modification at two-loop 
order where new flavor combinations are possible. 



Eq. (18 



reduces to the matching for the single-hadron 
Hiq. ^ . Secondly, we verify explicitly at one loop 

is equal to that 
Pi. 



FJF in 

that the anomalous dimension of Q i 
of the jet function J,, by using Eq. (jlSh and Eq 



This is a straightforward calculation that also requires 
the one-loop anomalous dimension |18j 



/' 



d jg\E,R,u,ii) _ ,,, 



dfi 2(2tt) 2 



y? t >(E,R,n)8(l-u) 

—~PjiW » 



(22) 



where jj is the jet function anomalous dimension, as well 
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as the RG equation of dihadron FFs [101 IH1 135] 



with anomalous dimension 



(23) 



E 



du a s ([i) 



7T J V 



Z\ z 2 



u u 



^>J u(l-u) tt Ku'^J a ^)ll-w 



III. A FULL NONPERTURBATIVE 
JET CHARGE DISTRIBUTION 

In this section we will take a different approach. Our 
starting point will be a nonperturbative distribution 
Di(Q, k, /i) for the charge Q of a parton of flavor i for 
a given k. (This is not to be confused with the frag- 
mentation function D^, and can be distinguished by its 
arguments.) We assume that these distributions are nor- 
malized, 



J dQDi(Q,K, fi) = 1 



(24) 



In analogy to Sec. |TTJ we will calculate the RG evolution 
of Di{Q, k, /i) and the corrections from the jet algorithm, 
which we described by a Gi(E, R, Q, k, /i). 

We start by observing that the one-loop RG evolu- 
tion consists of splittings i — > ja(ij), with a(ij) given 
in Eq. (15). The charge is the sum of the charge of the 



branches, 



Q = z K Q 1 + (l- z) K Q 2 . 



(25) 



where the rescalings z K and (1 — z) K are necessary because 
momentum fractions in the branches are taken with re- 
spect to their initiating parton. This suggests the follow- 
ing structure for the renormalization, 

1 



Dr e (Q,K,fx) 



2S 

j 



JdQi dQ 2 dz (z, fj,) 



x Dj(Q l ,K,fx)D a ^(Q2,K,(i) 
x5[Q-z K Q 1 -{l-z) K Q 2 }. 
From the partonic one-loop calculation, we find 



Zij{z,fJ>) 



2S ij S(l-z) + ^--P ji (z) 



(26) 



(27) 



Eq. ( 26 ) requires regulating the splitting functions for 
z — > 0, which may be obtained from the familiar z — > 1 
regularizations, using P gg (z) = P gg (l — z) and P gq {z) = 
P qq (l — z). Taking the /x-derivative of Eq. ( |26| ), we find 
the following one-loop RG evolution of the charge distri- 
bution 

fi-^- Di(Q, n,fi) 

j 

x Dj(Qi, k, n)D a (ij)(Q 2 , k, fi) 
x6[Q-z K Q 1 -{l-z) K Q 2 ], (28) 



Pji{z). 



(29) 



A nontrivial property of Eq. ( 28 ) is that it preserves the 



normalization in Eq. (24). Taking the appropriate mo- 



ments, Eq. ( 28 ) reduces to the evolution for the average 
and width of the charge distribution in Eqs. ^ and ( 17 ). 
One advantage of the approach in this section is that it 
does not require multihadron FFs for higher integer mo- 
ments. It also allows us to describe non-integer moments 
of the jet charge distribution, for which there is no de- 
scription in terms of multihadron FFs. 



The generalization of Eq. ( 28 ) to n- loops is expected 
to be given by 

d 1 f r n+1 

^-7- Di(Q, k, fx) = / Y[dQ m dz m D jm (Q m , k, n) 



x 51 1- 



n+l 



n+1 



m—1 m— 1 



(30) 



This becomes increasing nonlinear, but the nonlinearities 
are of course loop suppressed. 

In analogy to the fragmenting jet function, we intro- 
duce Gi(E, R, Q, k, fj.). This is the jet charge distribution, 
which takes the jet restriction into account. Similar to 
the renormalization in Eq. ( 28 ) , wc find that the one- loop 
matching is given by 

Gi(E,R,Q,K,(j,) = 7^^2 JdQidQ 2 dzJij(E,R,z,n) 
j 

x D j (Q 1 , k, n)D a (ifi(Q 2 , k, n) 

x 5[Q - z K Q 1 - (1 - zfQ 2 ] . (31) 

The matching coefficients Jij are the same as for the 
FJF, but now also need to be regulated for z — > 0. This 
regularization may be obtained from the z — ¥ 1 regular- 
izations, using J gg {E,R,z, ji) = J r ffg (i?,i?, 1 — z, ju,) and 
J qg {E,R,z,[i) = J qq (E, R,l - z,fi). To "cancel" the 
factor of 1/2 in Eq. (31) at tree- level, we modify the 
tree-level matching coefficients at z = 0: 



J$(E,R,z,ii) = 5(l-z), 
J$(E,R,z,n) = 5(z), 
jW(E,R,z,n)=5(z)+5(l-z), 
J$(E,R,z,n)=0. 



(32) 



d/x 



- / dQi dQ 2 dzj®(z,[i) We have checked that Eq. (31) preserves the normaliza- 



tion in Eq. (24) 



dQGi(E,R,Q, k, n) = 1. 



(33) 
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Taking moments of Eq. (31 1, 

(Ql) = J dQQGi(E,R,Q,K,ft) , 

((Qk) 2 ) = J d QQ 2 Q l {E, R, Q, k, ft) , (34) 

results in equations that are consistent with the expres- 
sions for the average jet charge and the width of the jet 
charge distribution in Eqs. ^ and (13 1. 



IV. MONTE CARLO IMPLEMENTATION 



An interesting feature of the approach in Sec. |III| is 
that it can be exactly recast in terms of a Monte Carlo 
(MC). Specifically, the solution to the nonlinear differ- 
ential equation in Eq. (28) is naturally obtained using 
a parton shower. Refs. |44Tl46] worked on (improving) 
parton showers in the context of SCET, by studying the 
matching and RG evolution of the hard function, which 
describes the short distance physics. Our approach is or- 
thogonal to this: we are only attempting to describe jet 
charge (and similar observables) but at higher precision. 
Our MC consists of a (possible) perturbative splitting at 
the jet scale fi ~ 2£?tan(i?/2), a parton shower connect- 
ing the jet scale and the hadronization scale [to ~ 1 GeV, 
and a hadronization step at ftQ. We will only keep track of 
the parton type and momentum fractions in the shower, 
which are the variables relevant for jet charge, and thus 
do not attempt to give a fully differential description of 
the final state. 

Starting with a quark or gluon with energy E, the first 
step in generating the jet consists of (the possibility of) 
a perturbative splitting at the jet scale, corresponding to 
Eq. (31 1. The splittings are described by the following 



probability densities, 
P[q->q(z)g(l-z)] = 

P[9^ g(z)g(l-z)] = 

P[9 -> q(z)q(l-z)] = 



J qq (E,R,z,n) 



2{2-KfJ q {E,R,ft) 

Jgg(E,R,Z,fl) 

i(2nyj g (E,R,ft) 

Jgq(E, R,Z,fl) 

2(2tt)3 J g (E,R, ft) 



0<z< 1-5, 

5 <z< 1-5, 

<z< 1, 
(35) 

where ft ~ 2£'tan(i?/2) to avoid large logarithms. These 
weights are negative, as is common at NLO, so the phrase 
"probability" is not meant in the classical sense. The ad- 
ditional factor of 1/2 for g — > gg is due to identical parti- 
cles. Because of the singularities in the Jij (regulated by 
plus distributions) we introduced a cut-off 5. This can 
be thought of as a resolution parameter, and defines the 
no-splitting probability, 



P[q -> q] 



J qq (E,R,z,ft) 



{ _ s dZ 2(2tt)3 J q (E,R,pL) 



P[g 



»' = /' 

Ji-, 



dz 



JggjE, R,Z,ft) 

2{2vf J 9 {E,R,[i) 



(36) 



5 dlnp 


{QD VI 


0.2 0.8 
0.2 0.4 
0.1 0.4 
0.05 0.2 
0.05 0.05 


0.159 0.387 
0.152 0.386 
0.146 0.370 
0.146 0.359 
0.142 0.358 


analytic 


0.142 0.353 



TABLE I. The average and width of the jet charge distribu- 
tion obtained from the JetFrag MC described in this section, 
for various values of the parameters 5 and dln/i, compared 
to the analytic calculation. We use the Gaussian toy model 
described in the text for the nonperturbative input. 



For sufficiently small values of 5 the result becomes inde- 
pendent of 5, but smaller values also increases the com- 
putation time. 



We subsequently carry out the evolution in Eq. ( 28 ) be- 
tween the jet scale and the hadronization scale by solving 
this equation iteratively, with step-size d In ft 



Di(Q, n, fie dlnfJ ") — Di(Q, k, ft) 



(37) 



£ fdz^-Pj^dlnft 



x fdQ 1 dQ 2 5[Q-z K Q 1 -(l-z) K Q 2 } 
x Dj(Qi, k, ft)D a (ij)(Q 2 , k, ft) . 



This builds up a parton shower, where for each step d In fi 
the splitting probability densities are 

P[q^q(z)g(l-z)] = P qq {z) A ln/i < z < 1-5, 

IT 

P[g -> g(z)g(l-z)} = P ss (z)dln M 5 < z < 1-5, 



P[g q(z)q(l-z)} = ^ P qg (z)dlaft 



0< z< 1, 
(38) 



and the corresponding no-splitting probabilities are 
P[q ->q\ = l 

P[g^g] = l+ I dz Pg g (z) d In ft. ! :-!<)') 



dz s ^ P qq [z) d ln/i , 



'1-6 

1+ [ dz 
J is n 



7T 

a s (fi) 



Note that the resolution parameter 5 used here can in 



principle be different from the one in Eqs. ( 35 ) and ( 36 1 . 

The description up to this point has been fairly generic 
and does not rely on the observable. However, now we 
need to use that in the shower the charge is the z K - 
weighted sum of the charge of the branches, and that 
the charge distribution of the branches is sampled over, 
as described by the last two lines of Eq. (37). We will 



only sample over the charge distributions of the partons 
at the end of the shower, randomly assigning them a 
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JetFrag example for k T -like, E=100 GeV, R=0.5, k=\ 

1.5r 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
Hadronization model: 


1 1 1 1 1 1 1 1 1 | 


Gaussian 






- Step function 


J~" — H 




i 1 — 

i 

i 

i , — i 

i 
i 

LI 


i_ 


-i 
i 
i 

1 — | i 



1.5 



-0.5 



0. 



0.5 



1.5 



FIG. 4 

Monte Carlo described in Sec 



Jet charge distribution obtained using the JetFrag 
We use the Gaussian (or- 



IV 



ange solid) or step function (green dashed) toy model de- 
scribed in the text for the nonperturbative input. 



D g {Q, k, no) and D q (Q, tt, ho) with the same averages 
and widths. From Sec. |n] we know that the resulting 
jet charge distribution will have the same average and 
width. However, the shape of the distribution is quite 
different, as shown in Fig. [4j The parton shower does 
turn the step function into a more Gaussian distribution, 
which can presumably be understood as a consequence of 
the central limit theorem. To the extent that a Gaussian 
description of jet charge suffices, the analytical approach 
of Sec. [n] is of course much simpler. 

The procedure described in this section can directly be 
extended to observables similar to jet charge, by replac- 
ing Di{Q, k, no) with the appropriate function for that 
observable. There are no obvious obstacles in extending 
this approach to n-loop order, although this will involve 
I — > n + 1 splittings [see Eq. (30)]. It will be interesting 
to see if a similar approach can be employed for more 
general track-based jet observables that are also sensi- 
tive to soft radiation, such as the track mass of a jet. We 
leave this question for future work. 



charge using Di(Q, n, Ha) as a probability distribution. 
This is our hadronization "model" . It is perhaps sur- 
prising that such a simple approach to hadronization is 
possible, compared to the string fragmentation models 
used in Monte Carlo programs, but this is because we 
restrict ourselves to a specific observable. By weighting 
the charges of these final partons with their momentum 
fractions, we obtain the jet charge for this event. Gener- 
ating a sufficient number of events (around I0 5 to I0 6 for 
statistical errors at the percent level) yields a numerical 
calculation of the full jet charge distribution. 

As proof of concept, we have turned this description 
into a computer program, which we call the JetFrag 
Monte Carlo. This generates unweighted events, but has 
so far only been minimally optimized. For simplicity 
we restrict ourselves to k = I, and adopt a simple toy 
model for the nonperturbative input. Specifically, we as- 
sume that D g (Q, k, Ho) is Gaussian with average and 
width 0.6 at the scale Ho = 1 GeV. Similarly we take 
D q (Q,K,Ho) — Dq(—Q,K,Ho) equal for all fiv^j quark 
flavors, and describe it by a Gaussian with mean 0.2 and 
width 0.4. In Table|!J we show results for the average and 
the width of the jet charge distribution, obtained using 
this MC program for various values of the parameters 8 
and dln/x. We compare this to the analytic results, to 
get an idea of how small these parameters need to be for 
a reliable description. 

The histogram for the corresponding jet charge distri- 
bution is shown in Fig. [4] for 6 = d In /i = 0.05. Since 
we start with nonperturbative input that is Gaussian, it 
is not surprising that the resulting jet charge distribu- 
tion looks very Gaussian again. As an illustration we 
therefore also show the result using a step function for 



V. NUMERICAL RESULTS 

This section contains our numerical results for the av- 
erage and width of the jet charge distribution. We first 
study the perturbative convergence of our results, fol- 
lowed by a detailed comparison with Pythia. We con- 
clude by discussing the optimal choice for n. 



A. Perturbative Convergence 

We start by studying the perturbative convergence of 
our calculation. In Fig. [5] we show the results for the 
average charge of a A^-like quark jet with R = 0.5 and 
k = 1, as a function of the jet energy E. The curves are 

quark, like, R=0.5, k-\ 




E [GeV] 



5 We here ignore the b and c-quark thresholds. 



FIG. 5. The average charge (Qf) at LO and NLO for a fcr-like 
quark jet with R=0.5 and K = 1. The bands correspond to 
the perturbative uncertainties as explained in the text. 
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FIG. 7. The average charge for an anti-fcr quark jet is shown as function of the jet energy E for various values of re and R. The 
Pythia results for d (u) quarks are shown as squares (circles). The plots are normalized to 1 at E = 100 GeV and R = 0.5, 
which removes the dependence on the nonperturbative input and thus the quark flavor. 



normalized to 1 at E = 100 GeV, which removes the de- 
pendence on the nonperturbative parameter in Eq. ([6]). 
At LO we do not include the NLO jet algorithm cor- 
rections, i.e. we take Jij — 2(27r) 3 5fj . As Fig. [5] shows, 
the NLO corrections reduce the average jet charge by a 
non- negligible amount. 

The perturbative uncertainties are estimated by vary- 
ing the renormalization scale /i up and down by a factor 
of 2. To keep the normalization point fixed, we simulta- 
neously vary the scale in the normalization. We show un- 
certainty bands both with (darker) and without (lighter) 
this additional prescription in Fig. [5| In all the following 
plots we will use this additional prescription, which keeps 
the normalization point fixed and leads to smaller uncer- 
tainties. However, since these uncertainty bands do not 
quite overlap, they may be a bit too optimistic. In ad- 
dition, the prescription causes the NLO band to be only 
slightly narrower than the LO result. (Neither of these 
issues are present for the lighter uncertainty bands.) 



In Fig. 6] we study the convergence of {(Q l K ) 2 ) for i = 
q,<7, which enters in the width in Eq. (10). We can no 



longer completely remove the nonperturbative input by 
normalizing, because of the mixing between quarks and 
gluons. We therefore make an assumption for 



{{QlY) 
{{Ql) 2 ) 



at fiQ = 1 GeV, 



(40) 



which we for simplicity take equal for all five light quark 
flavors. The solid curves and uncertainty bands corre- 
spond to p = 1 and the dotted curves in Fig. [6] corre- 
spond to p = 2. We find again that the convergence is 
reasonable. The mixing causes the width to reduce more 
slowly as function of E. (For quarks the effect of the 
mixing is stronger if p is larger, whereas for gluons it is 
the opposite way around.) 
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K 


M-quark 
Pythia DSS AKK08 


d-quark 
Pythia DSS AKK08 


s-quark 
Pythia DSS AKK08 


0.5 
1 

2 


0.271 0.237 0.221 
0.144 0.122 0.134 
0.055 0.046 0.064 


-0.162 -0.184 -0.062 
-0.078 -0.088 -0.046 
-0.027 -0.030 -0.027 


-0.196 -0.504 -0.123 
-0.108 -0.214 -0.064 
-0.043 -0.064 -0.024 



TABLE II. Average charge of an e + e anti-fey jet with E - 
due to its poor charge separation. 



100 GeV and R = 0.5. We did not include HKNS in this comparison 



ft 


u-quark 
Pythia HKNS DSS AKK08 


d-quark 
Pythia HKNS DSS AKK08 


s-quark 

Pythia HKNS DSS AKK08 


gluon 

Pythia HKNS DSS AKK08 


0.5 
1 

2 


0.341 0.862 0.734 

0.242 0.383 0.333 0.373 

0.136 0.155 0.143 0.157 


0.338 0.785 0.707 

0.236 0.339 0.313 0.340 

0.127 0.134 0.131 0.129 


0.336 0.549 0.674 

0.237 0.225 0.314 0.298 

0.132 0.093 0.139 0.130 


0.356 0.813 0.773 

0.199 0.314 0.300 0.381 

0.073 0.117 0.109 0.127 



TABLE III. Width of the jet charge distribution Y\ for an e + e" anti-fc T jet with E = 100 GeV and R = 0.5. The dihadron 
contribution is calculated assuming no correlations, using Eq. (42 I. 



B. Comparison with Pythia 



2. Comparison using Fragmentation Functions 



1 . Setup 



In this section we compare our calculation for the 
average and width of the jet charge distribution with 
Pythia Q2] [T3]. Pythia results for quark jets are 
obtained from the process e + e~ — > 7/Z — > qq, where 
•/s = 2E leads to jets of roughly the desired jet energy. 
The gluon jets are obtained from pp — > gg by taking 
yjs = 2E and requiring that the outgoing gluons have a 
minimum px ~ E. We cluster the jets using the e + e~ 
version of the anti-fc^ algorithm in Fast Jet [47], and 
only keep jets that have at least 95% of the desired en- 
ergy 



quark and gluon, anti-& r , 7?=0.5, k=\ 




300 400 
E [GeV] 



600 



FIG. 8. The size of the quark/gluon mixing in the width of 
the jet charge distribution. 



First we compare the average jet charge, using the in- 
put from the HKNS g5], DSS gSJ |SD] and AKK08 [ST] 
fragmentation function sets collected in App. [5] with 
the result obtained from Pythia at E = 100 GeV and 
R = 0.5. As Table [TT] shows, Pythia is in reasonably 
good agreement with the FF results, given the large un- 
certainties from the FFs, for which we take the spread 
between the FF sets as an estimate. The main reason 
for this spread is poor charge separation, as discussed in 
App.[B} 

Since the knowledge of dihadron FFs is very limited, 
we need either input from data or further assumptions to 
make predictions for the width. For example, recently di- 
hadron FFs have been calculated in the N JL model [S5] . 
Here we will assume that the two hadrons are uncorre- 



lated. Since the evolution in Eq. (17) and the matching 



in Eq. ( 13 1 generate correlations, we need to be more spe- 
cific: we assume that the dihadron FFs are uncorrelated 
at the low scale /xq = 1 GeV, 



d'; iM ( Zi ,z 2: 

which leads to 



Mo J 



D^(zi,im,)D^{z3,Ho), (41) 



Df Q (K,K ) n ) = [Df(K,^ )} 2 . (42) 

From the tables with nonperturbative parameters in 
App.|Bj we then see that the dihadron contribution D®® 
is typically (much) small er th an the single hadron con- 
tribution D? 2 ■ In Table 



III 



we show the width of the 
jet charge in Pythia, as well as the width obtained us- 
ing Eq. (42). For k — 0.5 and 1 the width obtained 



from Pythia is consistently smaller than that obtained 
from FFs. Of course this could indicate significant di- 
hadron correlation effects and not necessarily a problem 
with Pythia. As the knowledge of (dihadron) FFs pro- 
gresses, the above analysis could start putting constraints 
on parameters (or tunes) of Pythia. 



11 




FIG. 9. The width of the anti-fcr quark jet charge distribution as a function of the jet energy E for various values of re and 
R. The Pythia results for d (it) quarks are shown as squares (circles). The plots are normalized to 1 at E = 100 GeV and 
R — 0.5, since they are almost independent of the quark flavor. 




FIG. 10. The width of the anti-fcr gluon jet charge distribution as a function of the jet energy E for various values of re and 
R. The Pythia results are shown as circles, and the plots are normalized to 1 at E — 100 GeV and R = 0.5. 



3. Perturbative Comparison 

Next we compare the dependence on E and R in 
Pythia with our calculation, for which we take Pythia's 
values at E = 100 GeV and R = 0.5 as input. In the left 
panel of Fig. [7J we compare the average jet charge as 
a function of the jet energy E for u and d quark jets 
and k = 0.5, 1, 2. Since our calculation predicts the same 
shape for u and d quark jets, we normalize the plot. This 
also removes the dependence on the non-perturbative in- 
put describing the hadronizatiorj^] The Pythia result 
has the same features as our result and agrees fairly well 
for k = 0.5 and 1. For larger values of k, Pythia starts 
to yield a distinctly lower average jet charge. In the 



We observe this in Pythia as well: the average jet charge before 
and after hadronization is the same up to an overall factor. 



right panel of Fig. [7J we compare different values of R 
for re = 1. Pythia predicts a slightly lower average jet 
charge but the agreement is quite reasonable. 

We have investigated the size of the nonlinearities and 
mixing in the evolution and fixed-order corrections to the 
width of the jet charge in Eqs. (13) and (17). The effect 
of the nonlinearities is less than 1% until you get down to 
energies of only a few GeV, so they are irrelevant. How- 
ever, the mixing effect is quite significant. In Fig. [8] we 
show the (normalized) widths of the jet charge for quark 
and gluon jets, with and without the mixing contribu- 
tion. 

In Fig. [9j we study the width of the quark jet charge 
distribution as a function of the jet energy. We again 
show normalized results, since these are almost the same 
for the different quark flavors. In the left panel we com- 
pare our calculation to Pythia for several values of re. 
The rise in our prediction for re = 0.5 is due to the 
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large mixing contribution from gluons (which was not 
included in Ref. [H]). Pythia's decrease for k = 0.5 in- 
volves a cancellation between the parton shower and the 
hadronization model, since before hadronization Pythia 
also predicts a (relative) increase. At the same time the 
C(A 2k ) corrections due to soft radiation are the largest 
for small values of k. (Since soft radiation does not affect 
the average jet charge, this would not be in contradiction 
with the agreement seen in Fig. [7]) At this point it is not 
clear if the discrepancy indicates a problem with Pythia 
or our calculation. For n = 1, 2 the difference between 
our calculation and Pythia is smaller, though Pythia 
still predicts a distinctly lower width. In the right panel 
we perform the comparison for different values of R with 
k = 1. The dependence on R is very similar, but Pythia 
is again slightly lower. 



In Fig. 10 we study the width of the gluon jet charge. 
There is reasonable agreement for k = 1 and 2. For 
k = 0.5, we find that the gluon width barely changes 
due to the mixing with quarks, whereas Pythia yields a 
distinctly smaller width. As the right panel shows, the R 
dependence in Pythia exhibits the same general features 
as our result, but the effect is smaller. These points also 
suggest that the gluon jets in Pythia do not depend on 
E and R solely through the combination 2i?tan(i?/2) [as 
predicted by our calculation and observed for quark jets 
in Fig. |. " 

Based on these comparisons it seems reasonable to use 
Pythia for first jet charge studies. Precision studies will 
presumably lead to discrepancies with Pythia, though 
this can of course be alleviated by retuning. By con- 
trast, our higher-order calculation can be systematically 
improved to match the experimental precision. 



C. Optimal Choice for n 

Theoretically, the optimal choice the value for 

k where the peak of the jet charge distribution is best 
separated from zero. More precisely, K,^ IS where 



{Ql? 



(43) 



attains its maximum. Experimental considerations will 
of course affect this choice, since e.g. smaller values of k 
increase the (unwanted) sensitivity to soft radiation from 
the initial state or other jets |llj . Though the maximiza- 
tion of 77 is a nonperturbative question, we can study 



E(GeV) 


u d s 


25 
100 
400 


0.29 0.29 0.32 
0.28 0.25 0.31 
0.26 0.25 0.29 



TABLE IV. Theoretical optimal choice k» obtained from 
Pythia, for e + e~ anti-fcr jets with R = 0.5 and the indi- 
cated flavor and jet energy E. 




FIG. 11. Factor that controls the speed with which k* reduces 
as /1 ~ 2_Etan(i?/2) increases. Below k ~ 0.55 the reduction 
slows down. 



how depends on /1 ~ 2i?tan(i?/2). Neglecting NLO 
jet algorithm corrections and mixing with gluons, 



d?7 a s (/j,) 



dln/^ 
This leads to 
d/«* 



■2P„{k) - P qq (2 K )} V . 



(44) 



d In /i 



7T 77" (K* 



[P'M-PU^*)}, (45) 



so as we increase E and i?, K * IS reduced. The factor 
containing the splitting functions is plotted in Fig. 11 
and, interestingljj^J has a minimum around n ~ 0.55. 
Once IS below this value, its scale dependence will 
be reduced. In addition, a s (/z) reduces as /i increases. 
In table |IV| we show some results for obtained from 
Pythia. Since the values of sis &rc all well below 0.55, 
it is not surprising that there is little dependence on the 
jet energy. It should also be noted that 77 is fairly flat 
in the vicinity of k 4 , so a somewhat different value of 
may work almost as well. 



VI. CONCLUSIONS 

We have presented in detail the calculation of the jet 
charge distribution. This takes a particularly simple form 
for the average and width of the jet charge distribution, 
which are the experimentally most relevant parameters. 
The nonperturbative coefficients that enter in the aver- 
age and the width are related to moments of (dihadron) 
fragmentation functions. Since these currently still have 
large uncertainties the agreement with Pythia is rea- 
sonably good. Pythia tends to predict a width that 



7 By contrast, P' gg (K,) — P' gg {2K) ~ —1/k 2 , suggesting that any z- 
weighted property of gluon jets is quite diluted at high energies. 



13 



is smaller than those obtained from fragmentation func- 
tions when dihadron correlations are neglected, but this 
could of course be due to such correlations. We also 
compared our perturbative calculation with the shower- 
ing in Pythia, which agree reasonably, suggesting that 
Pythia suffices for initial studies of jet charge. The re- 
sults in this paper can be systematically improved by 
including higher order corrections, power corrections or 
updated (dihadron) FFs. 

There are various choices that enter in the jet charge, 
such as the jet algorithm, R and the weighting-power k. 
Jets with a smaller R retain a better jet charge signal, 
but also have an increased dependence on the jet algo- 
rithm. We studied the optimal choice k* of k for quark 
jets and found that it reduces as the jet energy increases. 
However, this energy dependence slows down for k* be- 
low 0.55, indicating that observables like jet charge can 
remain useful at high energies. 

We have also shown that our general calculation of jet 
charge can naturally be performed using a Monte Carlo- 
style approach. It is interesting that this is possible, given 
that standard Monte Carlo parton showers are limited to 
leading logarithmic order. We leave it to future work to 
investigate whether this can be extended to more general 
track-based jet observables that are also sensitive to soft 
radiation. 

In Ref. [TT] several potential applications of jet charge 
at the LHC were discussed. Here we confirmed through a 
detailed calculation that jet charge is theoretically under 
control. We therefore recommend a study of jet charge 
with LHC data as the natural next step. 
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Appendix A: Perturbative Coefficients 



1. One-Loop Splitting Functions 



where f3 = (11CU — 4n/7>)/3, is the lowest order coef- 
ficient of the QCD /3-function. In moment space [54] 



Pqq( K 

Pgg( K 

Pqg( K 



c f 



= c 



-2H(k + 2) + 

k 2 + 3k + 4 
k 3 + 3k 2 + 2k ' 

-C A [- 2H{k + 3) + 

k 2 + 3k + 4 

: If 



K+l K+2 



4k 2 + 6k + 4 
k 3 + 3k 2 + 2k 



k 3 + 6k 2 + 11k + 6 ' 

2r(K + i)r(K + 3) 

T(2k + 4) 



(A2) 



where H is the harmonic number function. The nonlinear 
contribution P qg decreases exponentially for large k, 



P gg( K ) = t f 



-(21ii2)k 



4^ 



1 + 



(A3) 



2. NLO corrections for kr-type jets 

The NLO matching coefficients Jij for the e + e~ ver- 
sion of fc^-like jet algorithms are given below, and we will 
also discuss its straightforward extension to pp collisions. 
At one-loop, where you have at most two partons, there 
is no distinction between the various e + e~ fc^-like jet al- 
gorithms. The jet restriction is simply 8 < R, where 9 is 
the angle between the two partons. This translates into 



s < 4z(l - z)E 2 tan 2 (E/2) 



(A4) 



where s is the invariant mass of the jet. For A^-like algo- 
rithms for pp collisions the corresponding jet restriction 
is y/ (Arj) 2 + (A0) 2 < R, where Arj and A(j> are the dif- 
ference in (pseudo)rapidity and azimuthal angle between 
the two partons. This is a more complicated restriction, 
as it is not rotationally symmetric around the jet axis. 
However, for narrow jets (R <C 1) it simplifies to [SS] 



s < 4z(l - z)p|tan 2 (i?/2) 



(A5) 



results can thus directly be extended to pp 
collisions by simply replacing the jet energy E by the 
transverse momentum px of the jet. 



The one-loop splitting functions are 
-1 + z 2 



P qq (z) = C F (^) + = C f 



n + z 2 a 

j r + -£(1- z) 

(1-z), 2 



P gq {z) = C F 
P ag (z) = 2C A 



z 

z 



(1-z), z J 2 



P qg (z)=T F [z 2 + (l-z) 2 }, 



(Al) 



K 


Jqq 


Jqg 




Jqq 


3gq 


0.5 


1.07 


-8.49 


-6.83 


-0.81 


-0.39 


1 


2.10 


-2.56 


0.02 


-0.64 


-0.14 


2 


3.65 


-1.01 


3.01 


-0.50 


-0.023 


3 


4.85 


-0.67 


4.51 


-0.43 


-0.0045 


4 


5.86 


-0.53 


5.64 


-0.38 


-0.00094 



TABLE V. Numerical results for the moments of the one-loop 
Ji^ (does not include the overall a s /iv and color factor). 
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K 


HKNS NLO 

D$ D'j D? 


DSS NLO 

5? D? 


AKK08 
D$ D Q d 5? 


0.5 
1 
2 


1.207 -0.807 -0.073 
0.420 -0.279 -0.062 
0.135 -0.089 -0.039 


0.302 -0.235 -0.642 
0.184 -0.132 -0.323 
0.087 -0.057 -0.121 


0.279 -0.079 -0.156 
0.199 -0.068 -0.095 
0.120 -0.051 -0.045 



TABLE VI. Nonperturbative parameters D®(k, fi = 1 GeV) for the average charge of a quark jet. 



K 


HKNS NLO 
5f 5f Df Df 


DSS NLO 

Df Df Df Df 


AKK08 
Df Df Df Df 


1 

2 
4 


0.676 0.540 0.113 0.680 
0.206 0.161 0.051 0.291 
0.050 0.038 0.016 0.108 


0.498 0.442 0.347 0.620 
0.165 0.143 0.130 0.260 
0.045 0.038 0.040 0.090 


0.188 0.150 0.099 0.447 
0.050 0.034 0.035 0.120 





HKNS NLO 


DSS NLO 


AKK08 


K 


Df Df 


Df Df 


Df Df 


1 


0.510 0.607 


0.618 0.639 


- 0.459 


2 


0.108 0.071 


0.145 0.081 


0.086 0.060 


4 


0.013 0.004 


0.022 0.008 


0.011 0.004 



TABLE VII. Nonperturbative parameters D^ (re, fi) that contribute to the width of jet charge distributions. For i = u, d, s, g 
the scale fi = 1 GeV and for i = c, 6 it is fj, = m c ^. For AKK08 most of the k = 1 moments are divergent, as denoted by 



Using the bare results in Ref. [T5J, we readily obtain, 
in the MS scheme, 



2(2tt) 3 



Jqg{E,R,z,n) 
2(2tt) 3 

2(2tt) 3 



PM-T.CpSil-z) L+C F J^(z) 



3 
-( 

2 

?M[P gg ( z )L + C F j£{z)], 
5 (l- z ) + ^El{c A L 2 8{l-z) 



P 3 g(z)-~PoS(l-z) L+C A J^(z) 



J gq (E,R,z,n) a a {n) 



2(2tt) 3 

Here L = ln[2£tan(i?/2)/^] and 



[P qg (z)L + T F J^(z)]. (A6) 



ln(l - z) 



+ 



- z /+ 
1 + z 2 , I 



+ (1 - z)ln(l - z) 

^2 



\nz + -(I - z) 
l-z 2 y ' 



l + (l-z) 2 



ln[z(l - z)] + 



5(1 -z) 



■^)=K 1 4^>v + 2(1+ '! )(1 -" b„ -:, 



z /+ 
2(1 -z + z 2 ) 2 

.2 i / 1 _\2i 



In ; 



24 



5(1 -z), 



J£{*) = W + (1 - z) 2 ] M*(l ~ ^)] + *(1 - *) • (A7) 

Anti-quarks have the same coefficients as quarks, and J qq 
and J qq i only start at two-loop order. We have checked 
these results using the sum rules in Ref. 18J. Since we 
evaluate the matching coefficients at fi ~ 2i?tan(i?/2), 
where L ~ 0, the one-loop contribution essentially comes 
from jfj 1 ■ We give numerical values for its moments in 



Table [v] Note that the nonlinear contribution J gq be- 
comes negligibly small for large k. 

The jet functions for fc^-like jets are [2"T] 



J q (E,R,n) = l 
J g (E,R,fj,) = l 



ol s (h)C. 



7T V 



o 3 r 13 3tt 2 

l — L^ 

2 4 8 



ot a (n) 



7T 

5 3vr 
24 ~ ~8 



M 2 - |L- 



23 
24 



(A8) 



Appendix B: Nonperturbative Coefficients from 
Fragmentation Functions 

In Table |VI| we show the nonperturbative parameters 
D9(k, fJ, = 1 GeV) for the average charge of a quark jet, 
using the HKNS 08], DSS 031501 and AKK08 [5T] frag- 
mentation function sets at NLO. The large differences be- 
tween the various FF sets are mainly due to poor charge 
separation. This is because we need the charge-separated 



combination D q — D q 



D h q 



Dg, whereas a lot of the 



data is e + e — > hX which only gives access to D q + D^. 
In particular, HKNS only uses e + e~ data in their analy- 
sis, so their quark/anti-quark separation relies crucially 
on assumptions. The large difference for the s-quark 
in DSS compared to the other FF sets is due to semi- 
inclusive DIS data that only they include [4TJ] , 

The results for , which is the contribution from 
single hadron FFs to the width of the jet charge distri- 
bution, is shown in Table |VII| Several of the k = 1 mo- 
ments for AKK08 are divergent (at z — 0) and denoted 
by a "-" . Since D® essentially depends on the combina- 
tion + D^ , the agreement between the different sets is 
much better. The gluon FFs are not as well known as the 
quark FFs, as is clear from the differences between the 
FF sets. We have also included the heavy quark flavors 
because they contribute through RG mixing. 
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